July 16, 2009 19:12 WSPC/INSTRUCTION 

KuznMikhOkrIJMPA.tex 



FILE 



International Journal of Modern Physics A 
© World Scientific Publishing Company 



REEXAMINATION OF A BOUND 
ON THE DIRAC NEUTRINO MAGNETIC MOMENT 
FROM THE SUPERNOVA NEUTRINO LUMINOSITY 



A.V. KUZNETSOV, N.V. MIKHEEV and A. A. OKRUGIN 

Division of Theoretical Physics, Department of Physics, 
Yaroslavl State University, Sovietskaya H, 
150000 Yaroslavl, Russian Federation 
avkuzn@uniyar.ac.ru, mikheev@uniyar.ac.ru, okrugin@uniyar.ac.ru 

Received (Day Month Year) 
Revised (Day Month Year) 

The neutrino helicity-flip process under the conditions of the supernova core is reinves- 
tigated. Instead of the uniform ball model for the SN core used in previous analyses, 
realistic models for radial distributions and time evolution of physical parameters in the 
SN core are considered. A new upper bound on the Dirac neutrino magnetic moment is 
obtained from the limit on the supernova core luminosity for vr emission. 
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1. Introduction 

Nonvanishing neutrino magnetic moment leads to the helicity-flip process where the 
left-handed neutrinos produced in the stellar interior could convert into the right- 
handed neutrinos being sterile with respect to the weak interaction, and this can 
be important e.g. for the stellar energy-loss. 

A considerable interest to the neutrino magnetic moment arised after the great 
event of the S'7V1987^4, in connection with the modelling of a supernova explosion, 
where a gigantic neutrino outflow defines in fact the process energetics. It means 
that such a microscopic neutrino characteristic, as the neutrino magnetic moment, 
would have a crucial influence on macroscopic properties of these astrophysical 
events. Too large outflow of right-handed neutrinos, produced due to the magnetic 
moment interaction, from the core would leave no enough energy to explain the ob- 
served neutrino luminosity of the supernova. Thus, the upper bound on the neutrino 
magnetic moment can be established. 

The neutrino helicity flip vi, — * vr under physical conditions corresponding 
to the central region of a supernova has been studied in a number of works (see, 
e.g., Refs. HHSl a more extended reference list is given in Ref. 0]). The process is 
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possible due to the interaction of the Dirac-neutrino magnetic moment with a virtual 
plasmon, which can be both generated and absorbed: 



In Ref.[TJ the neutrino helicity flip was described in terms of scattering by plasma 



core immediately after the collapse. However, the important polarization effects of 
the plasma on the photon propagator were not considered in that work. Instead, 
the photon dispersion was taken into account phenomenologically by introducing 
the so-called thermal mass of a photon into the propagator. The above-mentioned 
effects were considered more consistently in Refs. [2j |3j where the effect of a high- 
density astrophysical plasma on the photon propagator was taken into account 
using the thermal field-theory formalism. However, an analysis of works showed 
that they concerned only the electron component of the plasma, namely, only the 
channel v^e" — * v^er , and only the electron contribution to the photon propagator, 
whereas the proton component of the plasma was not considered at all. This seemed 
to be even stranger because the plasma-electron and proton contributions to the 
neutrino spin flip were estimated earlier ^ to be of the same order. 

A detailed analysis of the processes ((T|) , with neutrino- helicity conversion due to 
the interaction with both plasma electrons and protons via a virtual plasmon and 
with taking account of polarization effects of the plasma on the photon propagator 
was performed in Ref. [4] In particular, according to the numerical analysis, the 
contribution of the proton component of the plasma was not only significant, but 
even dominant. 

However, all those analyses HHH we re based on a very simplified model of the 
supernova core as the uniform ball with some averaged values of physical parameters. 
Moreover, the parameter values look, in modern views, rather too high than typical. 
It should be mentioned also that the improvement of the bound of Refs. [51 [3] with 
respect to the bound of Ref. Q] was based in part on the enlargement by the factor 
of 2 of the supernova core volume, while the core density was taken to be the same, 
p c ~ 8 x 10 14 g cm" 3 . This means that the core mass appeared to be in Refs. [51 [3] 
of the order of 3M , which is nearly twice as large as the mass of the supernova 
remnant believed to be typical. 

The aim of this paper is to make the estimation of the Dirac neutrino magnetic 
moment from the limit on the supernova core luminosity for vr emission by a 
more consistent way, taking some radial distributions and time evolution of physical 
parameters from some realistic models of the supernova core. 

For completeness, we consider here a general case of the magnetic moment matrix 
\i ViVj = fJiij (i.e. both diagonal and transition magnetic moments), where 2/j, i/j are 
the neutrino mass eigenstates. The neutrino states i>i with definite flavors i created 
in weak processes are the superpositions of the neutrino mass eigenstates: 



vl^vr + 7*, vl+7*-*vr- 



(1) 



electrons and protons [y^e 



VRe and v^p — > vrp, respectively) in a supernova 




(2) 
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where Uu is the unitary leptonic mixing matrix by Pontecorvo-Maki-Nakagawa- 
Sakata. It means that the value of the magnetic moment squared in all equations of 
Ref. |4] should be considered as an effective value. For the processes with the initial 
electron neutrino one should replace 



and similarly for the muon and tau initial neutrinos. 

The paper is organized as follows. In Sec. we give a clear illustration of the 
fact that neutrino scattering by protons dominates over their scattering by plasma 
electrons, basing on an analysis of a simplified case of the completely degenerate 
plasma, T = 0. Sec. [3] contains a summary of the procedure SI of obtaining the 
upper bound on the electron-neutrino magnetic moment from the «SW1987v4 data, 
in the frame of the uniform ball model for the supernova core. In Sec. 2] we make the 
estimation by a more reliable way, with taking account of radial distributions and 
time evolution of physical parameters, from realistic models of the SN core. The 
upper bounds are obtained on the combination of the effective magnetic moments 
of the electron, muon and tau neutrinos from the condition of not-spoiling the 
subsequent Kelvin — Helmholtz stage of the supernova explosion by emission of right- 
handed neutrinos during a few seconds after the collapse. 

2. Illustration: completely degenerate plasma at T — 

The comparison of the typical parameters of the supernova core, where the temper- 
ature is believed to be of order T ~ 15-30 MeV, while the electron and neutrino 
chemical potentials are rj e ~ 200-250 MeV and r\ Vc ~ 100 MeV, respectively, shows 
that the temperature is the smallest physical parameter. Thus, the limiting case 
of the completely degenerate plasma, T = 0, seems to give a reasonable estimate. It 
is remarkable that for the zero temperature limit the contributions from neutrino 
scattering by protons and electrons to the neutrino creation probability can be eval- 
uated analytically using Eqs. (20) and (21) and the corresponding formulas from 
Appendix A of Ref. [U 

It is appropriate to analyse the function T VR (E) defining the energy spectrum of 
right-handed neutrinos. In other words, this function specifies the number of right- 
handed neutrinos emitted per 1 MeV of the neutrino energy spectrum per unit time 
from unit volume of the central region of a supernova: 



a Hereafter we consider neutrinos as a quasiequilibrium gas described by the distribution functions 
fu(T, T)i/ e ) for the electron neutrinos and fu(T, 0) for the muon and tau neutrinos. This is believed 
to be a rather good approximation inside the SN core during a few seconds after the collapse. 



2 




(3) 



dn l/R 
dE 



(4) 
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The contribution of ultrarelativistic electrons to the function T UR (E) in the case 
T = can be obtained from the above-mentioned formulas of Ref. [4] in the simple 
form: 



9 9 

rS (E) = *±-fL {r, Ve ~ E) Ofo. - E) , (5) 

where E is the right-handed neutrino energy, fj, Ve is the effective electron neutrino 
magnetic moment ([3]) , = 2arj^/Tr is the squared mass of a transverse plasmon 
at T = 0, and 9{x) is the step function. 

The analytical expression describing the proton contribution turns out to be 
more complicated since it depends additionally on the proton mass. The plasma 
charge neutrality condition for T — takes the form n p — n e ~ and ensures that the 
electron and proton Fermi momenta are the same: k F e ^ = fcp . Then, the proton 
chemical potential coinciding with the Fermi energy is rj p = E F P ^ = ^Jm p + Ve 
and the proton contribution is expressed in terms of the proton Fermi velocity 
vp = k^/E^ = r] e /T]p — rj e / <J m p + r\\. As a result, the proton contribution can 
be expressed in the form: 

^E) = f^^±My), v--, o< y <i. (6) 

Here, the function <p p (y) has the different forms in two intervals: it is 

1 + -u F /3 

<Pp(y) = —, — —y> ( 7 ) 

1 — 1>F 



for < y < (1 — «f)/(1 + vf), and 



<Pp(y) 



1 



(1 - V F ) 



2 



1- \ 9 2 ' (l-y)(l + 2y) 
12 y z vp 



(8) 



for (1-Uf)/(1 + uf) <y < 1. 

Note that the formal turn to the limit m p — > 0, i.e. vp — > 1, in Eqs. ©-([5]) yields 
Vp(y) — * Ve(y) = (1 — y) #(1 — y), where the function ip e (y) can be introduced in 
Eq. ([5]) in complete analogy with Eq. ([6]). Thus, as expected, Eq. (|5|) for the electron 
contribution is reproduced. 

Figure 1 shows the plots of the function ip p (y) for vf = 1, Vf = 0.394, and vf 
= 0. The value v F = 0.394 corresponds to the effective proton mass m p ~ 700 MeV 
in a plasma with a nuclear density 3 x 10 14 g cm -3 (see Ref. El p. 152). The value 
vp = corresponds to the formal limit m p — > oo for which this function is also 
significantly simplified: (f p (y) -> (foc{y) =y6(l-y). 

The function T UR (E) defined in Eq. (j4|) determines as well the right-handed 
neutrino emissivity of a supernova core, i.e. the energy passed away by right-handed 
neutrinos per 1 MeV of the neutrino energy spectrum per unit time from unit 
volume: 
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y 

Fig. 1. Plots of the function <p p (y) for various up values. The dependence ip e (y) = (1 — y) for 
the electron contribution is reproduced for tip = 1 (dashed line). The value vp = 0.394 (solid 
curve) corresponds to the effective proton mass m p ~ 700 MeV. The case vp = (dotted line) 
corresponds to the limit of infinitely large proton mass. 



Q^= E ^ = ^n(E). (9) 

According to Eqs. ([4]) and ([9j), the right-handed neutrino emissivity is given by 
the formula 

2 2 4 
Lb Tfi 11 

Q»n = \J 3 ' v 3 iMv) + My)} ■ (io) 

The difference between the electron and proton contributions to the quantity 
given by Eq. (flO|) is illustrated in Fig. [2] It is clearly seen that the factor y 3 causes 
the increasing of the proton contribution to the emissivity. 



3. Uniform ball model for the SN core 

The spectral density of the supernova core luminosity via right-handed neutrinos is 
defined by the function T UR (E) as follows: 

^ = v ^ e = v ^ (E) = v y 3 (y, t) . (11) 

Here, V is the volume of the neutrino-emitting region, m 7 is the mass of a 



transverse plasmon, 



2a / , tt 2 T 2 



™? = 7(le 2 + ^|. (12) 
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Fig. 2. The function y 3 tp(y) defining the contributions from electrons (dashed line) and protons 
with m p ~ 700 MeV (solid line) and m p — > oo (dotted line) to the right-handed neutrino emissivity 
at T = 0. 




The function <^ nuro ) (y, T) introduced in Eq. (TT]) similarly to Eqs. © and flTUJ) can 
be extracted from Ref. |31 The function y 3 (^( num ) (y, T) is plotted in Fig. [3] for two 
values of the averaged temperature and for the electron and electron-neutrino chem- 
ical potentials rj e ~ 300 MeV and t\ Vf ~ 160 MeV. We neglected in our analysis SI 
the contributions of the processes with the initial muon and tau neutrinos. However, 
as will be shown below, these contributions appear to be essential. 

A comparison of Figs. [2] and [3] shows that taking of a nonzero temperature leads 
to a shift of the maximum of the energy distribution of the luminosity towards 
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higher energies of right-handed neutrinos. This additionally enhances the proton 
contribution. 

As a result, using the data on supernova SN1987A, a new astrophysical limit 
was imposed^ on the electron-neutrino magnetic moment: 



This is a factor of two better than the previous constraint!^' We have to remind, 
however, that both the previous and this improved bound on the electron-neutrino 
magnetic moment were based on a very simplified model of the supernova core as 
the uniform ball with some averaged values of physical parameters. In addition, the 
parameter values were set too high. For example, the upper limit 1.5 x 10 _12 ^b 
in Eq. (fl~3|) corresponds to the SN core temperature 30 MeV, while the limit 
0.7 x 10 _12 ^b corresponds to the temperature 60 MeV. As is seen from Fig. [3j 
the right-handed neutrino emissivity grows with temperature very rapidly. How- 
ever, according to recent simulations of the SN explosion, the temperature values 
inside the SN core are believed not to exceed 40 MeV, see e.g. Fig. 4. Anyway, 
taking account of the radial distribution of physical parameters inside the SN core 
would give more solid results. 

4. Models of the supernova core with radial distributions of 
physical parameters 

In this section we make the estimation of the upper bound on the Dirac neutrino 
magnetic moment by a more reliable way, with taking account of radial distributions 
and time dependences of physical parameters from realistic models of the SN core. 
Here we consider the models in the inverse chronology. 

4.1. The recent model of the O-Ne-Mg core collapse SN 

The very recent model was developed by H.-Th. Janka with collaborators who pre- 
sented us the results of their simulations^ of the O-Ne-Mg core collapse supernovae 
which were a continuation of model simulations of Refs. [3 [U The successful ex- 
plosion results for this case have rece ntly b een independently confirmed by the 
Arizona/Princeton SN modelling group EEO] w hich found very similar results. So 
we were provided with a model whose explosion behavior was comparatively well 
understood and generally accepted. 

We redefine Eq. (fTTj) . where, instead of multiplying by the volume of the 
neutrino-emitting region V, we integrate over this volume to obtain the spectral 
density of the energy luminosity of a supernova core via right-handed neutrinos: 



Here, taking the values defined in Eqs. (20) and (21) and the corresponding for- 
mulas from Appendix A of Ref. 131 we take account of their dependence on the radius 



/i„ < (0.7-1.5) x 10" 12 At B . 



(13) 




(14) 
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R [km] 

Fig. 4. The radial distribution for the temperature within the SN core at the moment t = 1.0 sec 
after the bounce, Ref. 6. 



R and time t. A comprehensive set of parameter distributions used in our estimation 
includes the profiles - of the density p, the temperature T, the electron fraction Y e , 
the fractions of electron neutrinos Y Ue , electron anti- neutrinos Yp c , and the fractions 
Y"„ x for one kind of heavy-lepton neutrino or antineutrino (y x — ^i, T ,^i,r), which 
are treated identically. The time evolution of the parameter distributions is calcu- 
lated^ within the interval until ~ 2 sec after the bounce. For the sake of illustration, 
we present in Figs. 4-6 the radial distributions within the SN core, from to 20 km, 
at the moment t = 1.0 sec after the bounce for the temperaturep' for the chemical 
potentials of electrons rj e and electron neutrinos r\ Vc (calculated on the base of the 
data of Ref. [6]), and for the proton nonrelativistic chemical potential r/* = rj p — m* N 
defining the degeneracy of protons (calculated on the base of the data of Ref. [6] and 
of the effective nucleon mass m* N in plasma, see Ref. p. 152). 

To analyse the influence of the right-handed neutrino emission on the SN energy 
loss, we also used the time evolution of the total luminosity of all species of left- 
handed neutrinos,^ presented in Fig. 7. 

Integrating Eq. ()14ll over the neutrino energy, one obtains the time evolution of 
the right-handed neutrino luminosity: 

oo 

W*)=2^ / dV J dEE 3 T VR (E). (15) 
o 

This is a novel cooling agent which would have to compete with the energy-loss via 
active neutrino species in order to affect the total cooling time scale significantly. 
Therefore, the observed SN1987A signal duration indicates that a novel energy-loss 
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R [km] 

Fig. 5. The radial distributions for the chemical potentials of electrons (solid line) and electron 
neutrinos (dashed line) within the SN core at the moment t = 1.0 sec after the bounce. 



100 



^> 50 

o 







ge-so 



-100 I — ■ — ■ — ■ — ■ — ' — ■ — ■ — ■ — ■ — ' — ■ — ■ — ■ — ■ — ' — ' — ■ — ■ — ' — lJ 
5 10 15 20 

R [km] 

Fig. 6. The radial distribution for the proton nonrelativistic chemical potential rj* — Tjp — 
within the SN core at the moment t = 1.0 sec after the bounce. 



via right-handed neutrinos is bounded by 

L UR < L UL . (16) 

Within the considered time interval until 2 sec after the bounce, one obtains from 
Eqs. (TT5)) . (fTB)) the time-dependent upper bound on the combination of the effective 
magnetic moments of the electron, muon and tau neutrinos. Assuming for simplicity 
that these effective magnetic moments are equal, one obtains the time evolution of 
the upper bound on some flavor-averaged neutrino magnetic moment JX V shown in 
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t [sec] 

Fig. 7. The time evolution of the total luminosity of all active neutrino species, Ref. 6. 




t [sec] 

Fig. 8. The time evolution of the upper bound on the neutrino magnetic moment within the 
time interval until 2 sec after the bounce (in assumption that the effective magnetic moments of 
electron, muon and tau neutrinos are equal). 



Fig. 8, where (i 12 = ^/(10 12 ^b)- 

As is seen from Fig. 8, the averaged upper bound tends to some value, providing 
the limit 

\± v < 2.4 x 10 12 /iB • (17) 
In a general case the combined limit on the effective magnetic moments of the 
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electron, muon and tau neutrinos is 

+0.71 (^ +M 2jj 1/2 <3J x 10-12 MB , (18) 

where the effective magnetic moments are dehned according to Eq. ([3]) . This limit is 
less stringent than the bound (fT3")) obtained in the frame of the uniform ball model 
for the SN core, but it is surely more reliable. Additionally, the upper bound on the 
effective magnetic moments of muon and tau neutrinos is established. 

4.2. Earlier models of the SN explosion 

The s imi lar procedure of evaluation was performed with using of the data of the 
modepJJ by R. Buras ct al. (2006) of the two-dimensional hydrodynamic core- 
collapse supernova simulation for a 15 M Q star. Namely, the radial distributions 
of parameters at the moments t = 0.2, 0.4, 0.6, 0.8 sec after the bounce in the model 
sl5Gio_32.a were taken from Fig. 40 of Ref . [TT1 Additionally, the fraction of electron 
neutrinos was evaluated as Y Vc ~ (1/5) Y e . Calculating the right-handed neutrino 
luminosity with those parameters and putting the limit (|16p . where the total lumi- 
nosity via active neutrino species L VL in that model can be taken from Fig. 42 of 
Ref. 1111 one obtains that the upper bound on the flavor-averaged neutrino magnetic 
moment p, v also varies in time as in the previous case. The time-averaged upper 
bound on \x v corresponding to the interval 0.4-0.8 sec, is: 

/2„ < 2.7 x 10- 12 /i B , (19) 

to be compared with the limit (fTTjl . 

Using the results of Ref. [12] by J. A. Pons et al. (1999) where the thermal and 
chemical evolution during the Kelvin-Hclmholtz phase of the birth of a neutron star 
was studied, taking the data from Figs. 9 and 14, we have obtained the time-averaged 
upper bound on /2„ for the time interval 1-10 sec of the post-bounce evolution in 
the form: 

p u < 1.2 x IO^Vb- (20) 

We also used the results of Ref.lHby W. Keil and H.-Th. Janka (1995) where the 
numerical simulations were performed of the neutrino-driven deleptonization and 
cooling of newly formed, hot, lepton-rich neutron star. Using the data presented in 
Figs. 3-9 on the SBH model (of the hot star with a "small" barionic mass), we have 
evaluated the time-averaged upper bound on ji v for the time interval 0.5-5 sec after 
the bounce in the form: 

fiu < 1.1 x 10- 12 /i B . (21) 

One can summarize that the upper bound on the flavor- and time-averaged neu- 
trino magnetic moment at the Kelvin-Hclmholtz phase of the supernova explosion 
occurs to be 

flu < (1.1-2.7) x 10- 12 /i B , (22) 
depending on the explosion model. 
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5. Conclusions 

The right-handed neutrino luminosity caused by the neutrino hclicity-flip process 
under the conditions of the supernova core, where the produced left-handed neu- 
trinos could convert due to the neutrino magnetic moment interaction into the 
right-handed neutrinos, being sterile with respect to the weak interaction, is reinves- 
tigated. Instead of the uniform ball model for the SN core used in previous analyses, 
realistic models for radial distributions and time evolution of physical parameters 
in the SN core are considered. The upper bounds on the flavor- and time-averaged 
magnetic moment of the Dirac type neutrino are obtained in those models, from 
the condition of not-affecting the total cooling time scale significantly: 



depending on the explosion model. 

In the recent paperpl the sterile right-handed neutrino luminosity was calcu- 
lated with taking account of radial distributions of the supernova core parameters, 
using the one-dimensional astrophysical code "Boorn"!!^ The supernova matter 
state parameters were calculated as the functions of coordinate and time during 
250 ms after the bounce. We should give a comment on Table 2 of Ref. [TU At first 
glance, the result for the right-handed neutrino luminosity, (0.5-1.1) xlO 50 erg/s, 
obtained in Ref. [14] is in agreement with the result (0.4-4) xlO 50 erg/s by R. Bar- 
bieri and R.N. MohapatraJ^ and contradicts to ouiSl result (3.8-22) xlO 50 erg/s. 
However, as is seen from the paper, the result of Ref. [14] for the minimal value of 
the right-handed neutrino luminosity, 0.5 xlO 50 erg/s corresponds to the effective 
temperature of the emitting matter T e ff ~ 10 MeV, while the minimal values of 
the luminosities of Refs. [1] and [4] correspond to the averaged temperature 30 MeV. 
Taking into account the strong dependence of the right-handed neutrino luminosity 
on the temperature, mentioned in Refs. [4] [TU one should conclude that the results 
of Refs. [T4l and l4l are in agreement, and both contradict to the result of Ref. 1. 
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ft, < (1.1-2.7) x KT 12 /i B , 



(23) 
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